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Abstract. In this paper we review the recent extraordinary progress in the 
development of a new quantum standard for resistance based on graphene. We 
discuss the unique properties of this material system relating to resistance metrology 
and discuss results of the recent highest-ever precision direct comparison of the Hall 
resistance between graphene and traditional GaAs. We mainly focus our review 
on graphene expitaxially grown on SiC, a system which so far resulted in the best 
results. We also briefly discuss progress in the two other graphene material systems, 
exfoliated graphene and chemical vapour deposition graphene, and make a critical 
comparison with SiC graphene in light of their possible application in quantum 
resistance metrology. 
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1. Introduction 

The discovery 8 years ago of the quantum Hall effect (QHE) in graphene sparked 
an immediate interest in the metrological community. The QHE is a fascinating 
macroscopic quantum effect occurring in two-dimensional conductors and relates the 
resistance quantum, h/e 2 only to the fundamental constants of nature, h, the Planck 
constant and, e, the elementary charge [H[2]- Although metrology successfully used the 
QHE for more than two decades to realise the resistance scale [3] , graphene is a material 
with properties like no other. Graphene, a single layer of carbon atoms in a hexagonal 
crystal lattice structure, is a truly two dimensional metal with a linear dispersion 
relationship characteristic for massless Dirac-type charge carriers jl]. The unique 
bandstructure of this semi-metal has both practical and fundamental implications. 
Firstly, the massless nature of the charge carriers leads to a Landau level spectrum 
with an energy gap between the first two levels which is around 5 times larger than that 
in semiconductor materials for magnetic fields around 10 Tesla. This implies that the 
QHE in graphene can be observed at much reduced magnetic fields and/or much higher 
temperatures [5]. Secondly, the marked difference in bandstructure and charge carrier 
characteristics between graphene and semiconductor systems allows for a demonstration 
of the universality of the quantum Hall effect through a rigorous test of the material 
independence of the value of Rk = h/e 2 , the von Klitzing constant. 

Theory predicts no corrections to the simple relation i?x = h/e 2 . The quantum 
Hall resistance is considered to be a topological invariant [6] , not altered by the electron- 
electron interaction, spin-orbit coupling, or hyperfine interaction with the nuclei. It has 
also been shown that the quantised Hall resistance is insensitive to much more subtle 
influences of the gravitational field [7]. Recently, a quantum electrodymical approach 
to charge carriers in a magnetic field has predicted a tiny correction to the von Klitzing 
constant of the order of 10~ 20 for practical magnetic field values [8]. However, the size 
of this predicted correction is about 8 to 10 orders of magnitude smaller than the most 
accurate measurement techniques available and therefore untestable. Nevertheless, the 
fundamental nature of the Hall resistance quantisation makes experimental tests of its 
universality of the utmost importance. Universality of h/e 2 will also strongly support 
the pending redefinition of the Si-units for kilogram and ampere in terms of h and e [9] . 

A comparison of the Hall resistance in two different substances does not prove the 
exactness of the relationship i?K = h/e 2 , however, material independence is a significant 
factor in establishing the fundamental nature of Rk- This material independence turns 
out to be rather difficult to establish. Indeed the characteristics of QHE samples must 
satisfy very stringent requirements pX] and in 30 years only silicon MOSFETs (metal- 
oxide-semiconductor field-effect transistors) and III-V's (GaAs/AlGaAs or InGaAs/InP 
heterostuctures) did so [3]. 

The first accurate measurements of the QHE in graphene were performed by 
Giesbers et al. [12] on exfoliated samples. The precision obtained in these measurements 
was at the part-per-million (ppm) level and limited by the high contact resistances 
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a) 




Figure 1. (a) Optical micrograph of a Hall bar used in the experiments, (b) Layout 
of a 7 x 7 mm 2 wafer with 20 Hall bars, (c) AFM image of Hall bar. Meandering lines 
are steps in the surface. (Adapted from Ref. [TP] ) 



together with a low maximum source-drain current which these samples could sustain 
before breakdown of the QHE occurs. A large measurement current determines the 
maximum signal-to-noise ratio and increasing this breakdown current is key to high- 
accuracy measurements. One established method of increasing the breakdown current 
is to increase the sample width [3] which is not easy to achieve with the exfoliation 
technique. 

A breakthrough came in 2009, when several groups within days of each other 
succeeded in growing large-areas wafers of epitaxial graphene by sublimation of SiC 
(SiC/G) with a quality good enough to observe quantised Hall resistance p3J [HI [151 
[T6| [T7] . In an indirect comparison with a GaAs/AlGaAs heterostructure device via an 
intermediate room-temperature standard resistor, Tzalenchuk et al. [15] were the first 
to demonstrate an accuracy of 3 parts in 10 9 (ppb) for the resistance quantisation in a 
large SiC/G device (see fig. [Q. The measurement system [18] was of identical design 
to that used by Giesbers et al. p2] and the key factor in the large improvement was 
the very low contact resistances which could be achieved and an order of magnitude 
increase in the breakdown current. The large breakdown current in SiC/G is not simply 
the result of a larger device but finds its origin in a charge exchange mechanism between 
the SiC substrate and graphene, leading to an unusually strong pinning of the quantum 
Hall state [191 ED]- This mechanism turns out to be stronger for lower carrier densities 
and by using a novel photochemical gating technique, Lara-Avila et al. [21] have been 
able to develop a graphene device with an extremely robust v — 2 quantum Hall state. 
Robustness, in this context means, invariance of the quantisation under changes of 
magnetic field, temperature or source-drain current and is a key requirement for any 
practical application. Recently, the same team undertook a direct comparison between 
SiC/G and GaAs, using an enhanced measurement system, demonstrating equivalence 
of i?H with a relative uncertainty of 8.7 parts in 10 11 [22]. 
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2. Quantum Hall effect in conventional 2DEGs 

When a conducting material is placed simultaneously in electric (E) and magnetic (B) 
fields, charge carriers (electrons or holes) moving with a velocity v experience a magnetic 
(Lorentz) force proportional to ±e \ E + v X By In the Hall bar geometry of fig. EH, 
carriers moving in the x direction will be deflected by the perpendicular magnetic field 
B z and accumulated at the edges of the sample. This accumulation of carriers gives 
rise to a potential build up in the direction transversal to carrier propagation (V xy ), the 
Hall voltage. Under these conditions the resistivity of the sample takes a tensor form, 
being different in the direction parallel or perpendicular to the direction of the current. 
We define the transversal resistivity p xy as: 

~^^xij ' Z / -i \ 

Pxy = — = — (1) 

^ XX ihG 

The constant of proportionality between the transversal resistivity and the applied 
magnetic field is called the Hall coefficient (Rh = 1/ne). The longitudinal resistivity of 
the sample p xx , on the other hand, is independent of magnetic field and can be obtained 
by measuring the voltage developed along the direction of the current and taking into 
account the geometry of the sample (length L and width W). 

wv xx 

Pxx = -r— (2) 

^xx 

In general, low-field Hall effect measurements are a simple yet powerful tool to 
study electronic properties of materials. The magnitude and sign of the Hall coefficient 
have been extensively used to characterize electronic materials in terms of carrier type 
(electrons or holes) and density n (per unit area). In addition, the Hall coefficient and 
the longitudinal resistivity can be combined to find the average carrier mobility in the 
sample as n — Rh/ p xy - 

The classic Hall effect picture breaks down when working with high mobility, two- 
dimensional systems at low temperatures and in strong magnetic fields [23] . Under these 
conditions, quantum effects are manifest as quantized steps in p xy , while simultaneously 
the longitudinal resistance vanishes, p xx = (fig. |2b). Steps and plateaux in the 
transversal resistance occur at universal values of p xy = h/ue 2 , with h the Planck 
constant, e the elementary charge and v and integer. This is the quantum Hall effect 
(QHE) [23]. 

The QHE can be understood as a consequence of Landau quantization [24] . Under 
the presence of a strong magnetic field, the size of the cyclotron orbit shrinks and 
becomes comparable to the wavelength of carriers. If the mobility in the system is high 

% For a two-dimensional conductor the 2D resistivity is denned as piu = p^d/I, with t the 
(infinitesimal) thickness of the sample. p2D has the same dimensions as the 3D resistance (il). In order 
to avoid confusion, the units of p2D are most of the times explicitly written as Cl/(L/W) = fi/square. 
Sometimes, however, the words resistance and resistivity for two-dimensional samples are used 
indistinctly. 
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Figure 2. (a) Magnetotransport measurements in a Hall bar geometry enable the 
electrical characterization of a material in terms of carrier type (electrons or holes), 
concentration and mobility, (b) Classic (quantum) Hall effect at low (high) magnetic 
fields in conventional, semiconductor-based, two-dimensional systems. 



enough so that carriers can complete a few cyclotron orbits before its momentum is 
relaxed, that is: 

uj~ 1 « r (3) 

then the size of the orbit can only take some allowed values, corresponding to an 
integer number times the electron wavelength. Recalling that the cyclotron frequency is 
given by u c = eB/m* and the carrier mobility can be expressed as [i = er/m, condition 
()3]) can be rewritten as B ^> yT 1 . Thus, quantum effects can in principle be observed 
with B = 1 Tesla if the mobility of the sample is around \i = 10, 000 cm 2 V" 1 s _1 . 

Quantized cyclotron motion modifies dramatically the density of states by breaking 
it into discrete levels, so-called Landau levels (LL), which are the allowed energies for 
cyclotron orbits under quantizing conditions (fig. [3^). Formally, LL are obtained by 
solving the Schrodinger equation for free electrons in a 2D systems, in the presence 
of electric and magnetic fields. The problem reduces to that of a harmonic oscillator 
shifted by the magnetic length £b = \fh/eB, with eigenvalues given by [251 Ell EE] : 

E N = hu c (N + 1) (4) 

with the reduced Planck constant h = h/2ir, the cyclotron frequency u c = eB/m, and 
iV an integer, zero included. In theory, each LL is a heavily degenerate delta function 
in the density of states filled with localized states undergoing cyclotron motion. In 
real samples, nevertheless, LL are broadened by the presence of disorder, (fig. [3^) due 
inelastic collisions (energy exchange) of localized electrons with impurities. 

The steps observed in p xy in QHE can be explained by considering that LL are 
filled with localized carriers undergoing cyclotron motion. As carriers undergo cyclotron 
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Figure 3. (a) Quantization of cyclotron motion leads to formation of Landau levels 
(LL) in the density of states (DOS). In theory the levels, spaced by fiu) c , are highly 
degenerate delta functions, but in real samples the presence of disorder broadens the 
levels, (b) In the absence of electrostatic gate the electron concentration remains 
constant and the quantum Hall plateaux are observed by increasing the magnetic field 
B, whenever the Fermi energy lies in between the the center of Landau levels. 



motion, they enclose a quantum of magnetic flux $o — h/e, and thus the amount of 
localized carriers per unit area till in the system can be found as : 

n LL = B/% = eB/h (5) 

With till known, the amount of completely filled Landau levels v (so-called filling 
factor) can be found by dividing the total density of electrons in the system, n, by the 
number of localized carriers, till. 

v = njn LL (6) 

When a LL is full the Fermi level lies in a gap between occupied levels and the 
filling factor v in Eqj6]must be an integer (fig. Ek). If we substitute the discrete density 
of states resulting from the formation of LL into the field- dependence of p xy for the 
classic Hall effect (Eq. CD), we find that p xy is quantized as: 

B B h 

Pxy = — = = — (7) 

en evn LL ue z 

A common interpretation for QHE is in the picture of extended states, carrying 
current without dissipation (zero resistance) along the edges of the sample, and localized 
states undergoing cyclotron motion in the bulk [27JI2E]- The origin of zero longitudinal 
resistance p xx is that extended states propagating in one direction of the sample are 
spatially separated from those carrying current in opposite direction, suppressing in this 
way backscattering. Maxima in p xx is observed every time the Fermi level crosses the 
center of a LL; plateaux in p xy and vanishing p xx are observed whenever the Fermi level 
lies in between the center of LL, pinned by localized states. Experimentally, quantum 
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Figure 4. Simplified schematic of the cryogenic current comparator bridge circuit. 
(Adapted from Ref. [18]) 



Hall plateaux can be observed by: a) fixing the magnetic field and varying the Fermi 
level of the sample with e.g. an electrostatic gate or b) fixing the Fermi level (fixed 
carrier concentration) and varying the magnetic field (fig. [3b). 



3. Resistance metrology 

3.1. Requirements 

In resistance metrology the value of the quantum Hall resistance, either at ~ 12.9 kQ for 
% = 2 or at ~ 6.4 kf2 for i = 4, needs to be compared with standard resistors at decade 
values typically in the range of 1 Q to 100 kQ. Room temperature measurement systems 
such as potentiometric or current transformer bridges can perform these measurements 
to a level of a few parts in 10 8 [3]. However, most top-level national measurement 
laboratories make precision measurements with a so-called cryogenic current comparator 
(CCC) bridge [29]. A CCC bridge can be used to determine the ratio of almost any pair 
of resistors with an accuracy of a few part per billion (ppb) [32]. If one of these resistors 
is a quantum Hall device, it can be used to realise the value of the other resistor in 
terms of the SI Ohm via the internationally agreed value of the von Klitzing constant, 
-Rk-90 [HQ] . If both resistors are quantum Hall devices, the CCC bridge can be employed 
to determine the universality of the quantum Hall effect |31j . 
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The bridge circuit for measuring the ratio of two resistors in terms of the current 
ratio determined by a CCC is schematically shown fig. HI Two isolated current 
sources "Master" and "Slave" separately drive current through resistors Rl and R2 
and associated windings Ni and N 2 on the CCC. The current ratio can be set via 
electronics to a few parts in 10 6 and this ratio is improved to a level of 1 part in 
10 11 by forming a negative feedback loop from the SQUID (Superconducting Quantum 
Interference Device) sensing the net flux in the CCC to the Slave current source. When 
balanced, the reading on the nanovoltmeter is exactly proportional to the resistance ratio 
(often a second servo loop is used to null the nanovoltmeter reading by injecting a current 
in a third, N3 small winding on the CCC [IE])- The main uncertainty components of 
the bridge are, the noise of the SQUID sensor, the Johnson noise of the resistors being 
measured and the current and voltage noise of the null detector. These components 
combine to give a typical measurement uncertainty of 1 ppb for a measurement time 
between 5 and 15 minutes |18j . 

In principle the value of the quantum Hall resistance is exactly e 2 /ih when the 
magnetic field is set to integer filling factor i. However, in practice there are many 
possible effects which can cause the Hall resistance to deviate from exact quantisation. 
One of the key factors is unwanted dissipation: if p xx is not exactly zero, the Hall device is 
not in the non-dissipative state and p xy will deviate from exact quantisation. Also if the 
contact resistances to the 2DEG are too large (typically > 10 Q), excessive dissipation 
will occur at the contacts which again can cause errors in R H . Over the last two decades 
the metrological community has evaluated all possible errors in the quantum Hall effect 
and composed a set of experimental guidelines which, when precisely followed, will result 
in a measurement uncertainty smaller than 1 ppb [U [H] . 

3.2. State-of-the-art 

When measuring two quantum Hall devices with a combined impedance of ~ 25 kfl 
the current-noise properties of the nanovoltmeter start to dominate the uncertainty 
budget. This problem can be alleviated by employing a second CCC in place of the 
nanovoltmeter (see fig. [5] where Rl and R2 are replaced by QHR samples SI and S2). 
The potential contacts on SI and S2 are closed in a loop via winding C on this second 
CCC (see fig. [5]). This device is configured with just a single winding to measure a 
current null rather than two windings to establish a current ratio. Janssen et al. [32] 
made a careful evaluation of all uncertainty components of this measurement system 
and demonstrated that a systematic uncertainty of only 1 part in 10 11 can be achieved 
for a direct 1:1 comparison of two quantum Hall devices. 

Figure [6] shows the Allan deviation for a 3.5 h measurement of graphene against 
GaAs at v = 2 of the CCC bridge in Ref. [32]. The Allan deviation is very convenient 
for a resistance ratio measurement involving periodic current reversal as it leads directly 
to the expected measurement resolution for a given measurement time. From the figure 
we see that the first data point for a ~ 40 s measurement gives an uncertainty of around 
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4.2 K Enclosure 



Figure 5. Schematic of the cryogenic current comparator bridge circuit used for high 
precision universality test. (Adapted from Ref. [22l 132) ) 

4 parts in 10 9 . The Allan deviation decreases as l/y/r expected for white noise. After 
3.5 h measurement time a relative uncertainty of 2 parts in 10 10 is achieved. There is still 
room for improvement here; if the two CCC's can be made to operate at their optimum 
noise performance, a 40 s measurements would give a relative uncertainty slightly better 
than 1 part in 10 9 and 6 parts in 10 11 after 3.5 h (blue dot and dashed line in fig. |6]). 

Even better performance for the universality test could be achieved by using four 
quantum Hall devices in a Wheatstone bridge configuration developed by Schopfer and 
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Figure 6. Typical Allan deviation of a 1:1 measurement of graphene against 
GaAs/AlGaAs heterostructure at v = 2 for a measurement current of 100 /xA. 
(Adapted from Ref. (32]) 



Poirier [33]. The QHR devices are connected in a series, parallel network using the 
multiple connection scheme developed by Delahaye and a CCC senses the imbalance in 
the bridge. A direct comparison of four GaAs devices has demonstrated an accuracy of 
4 parts in 10 11 and an uncertainty of a few part in 10 12 is conceivable [34] . 

3. 3. QH resistance in terms of the ohm 

The determination of the value of h/e 2 in terms of SI units is a rather complex 
experiment and not routinely performed. The derived SI unit for ohm is s _3 kgm 2 A~ 2 
and depends on four of the seven base units. In practice, the route to realise the ohm 
does not involve the ampere but rather relies on a clever electrostatics theorem by 
Thompson and Lampard [35]. The theorem states that the change in cross-capacitance 
of a specific four electrode device as a function of change in guard electrode only depends 
on frequency and length which both can be realised with very small uncertainty. From 
impedance we can get to resistance via a long and complex chain of a.c. bridges which 
at the moment can be done to a few parts in 10 8 |36j . 

The determination of Rk = h/e 2 in SI units also allows one to determine the value 

2 

of a, the fine structure constant, because a = In the SI, the permeability of 

vacuum fi and the speed of light c are fixed quantities with /i = 4n x 10~ 7 NA~ 2 
and c = 299792458 ms" 1 . Indeed, the first report on the quantum Hall effect was as a 
new and more accurate determination of a [23] . This determination of a can compared 
with other determinations such as high precision measurements and calculations of the 
anomalous magnetic moment of the electron, gyromagnetic ratio of protons or the mass 
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of neutrons which results in a combined value of the fine structure constant with an 
uncertainty of 3.2 parts in 10 10 . 

Early on after the discovery of the quantum Hall effect it was realised that the 
stability and reproducibility of the QHE was much better than the ability to realise its 
value in SI units. Therefore in 1990 it was decided to assign a constant value to Rk 
called -Rr-90 with value 25812.807 Q exactly without uncertainty [37]. This has allowed 
metrology laboratories to calibrate standard resistors with uncertainty typically in the 
range of 5 ppb or better [3]. 

4. Universality tests in conventional 2DEGs 

The fact that Rk = h/e 2 means that the quantum Hall effect should be universal, by 
which we mean that the value of Rk is independent of the sample properties (carrier 
density, mobility, dimensions, Landau level index and host material of the 2DEG). 
However, there is no quantitative theoretical model which proves this and the problem 
has essentially been addressed experimentally. The method described in the previous 
section is limited in accuracy by the complexities of an SI realisation [38] and a much 
higher precision can be achieved by directly comparing two quantum Hall devices against 
each other in a unitless null measurement. 

The first universality test was reported by Hartland et al. [31] who measured the 
ratio of the % = 2 plateau in a GaAs/AlGaAs sample and the % = 4 plateau in a Si- 
MOSFET sample and obtained R H (i = 2; G&As)/R H {i = 4; Si) = 2[1 - 0.22(3.5) x 
10 -10 ]. Later Jeckelmann et al. [39] improved on this result and reduced the uncertainty 
to 2.3 x 10~ 10 . Subsequently, the universality has been tested by varying the mobility, 
channel width and plateau index number [40]. From the results one can conclude that 
in semiconductor systems the universality has been established within a few parts in 
10 10 . 

5. Graphene, bandstructure, electronic properties and the quantum Hall 
effect 

Graphene is the first two-dimensional crystal available for experiments [U [2j HTJ 
142] . When compared to quasi-2D systems fabricated at interfaces (e.g. Si/Si0 2 
or AlGaAs/GaAs), graphene displays a peculiar linear E(/c) dispersion (instead of 
parabolic) that modifies substantially the magnetotransport picture. The following 
sections present an overview of the electronic structure of graphene and its implications 
for magnetotransport. 

5.1. Crystal and electronic Structure 

Graphene owes its two-dimensional nature to sp 2 hybridization of carbon atoms. The 
crystal backbone is formed by in-plane localized a bonds, between carbon atoms, and 
out-of-plane it electrons delocalized over the entire crystal. Each carbon atom has three 
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Figure 7. (a) Graphenc honeycomb structure (direct space); the unit cell contains 
two carbon atoms separated by clq = 1.42 A. (b) Reciprocal lattice showing the first 
Brillouin zone (FBZ); six points at the corner of the FBZ fall into two groups of 
equivalent points, denoted as K and K'. (c) Calculation of graphene bandstructure 
using a hoping parameter t = —3 cV; the bands are shown along the main 
crystallographic directions, including the K and K' points. 

nearest neighbors, separated by ~ 120° giving graphene its characteristic honeycomb 
structure. 

The graphene lattice is described by a unit cell that comprises two atoms, A and 
B, periodically arranged in a triangular lattice (fig. [7k). In real space, the primitive 
vectors are given by a\ = xa + yb and a 2 = xa — yb, a = 3a /2, b = y/3a /2 and the 
distance between nearest neighbors a = 1.42 A. The reciprocal lattice is constructed as 
K = MAi + NA 2 , where (M,N) are integers and the primitive vector in the reciprocal 
lattice Ai and A 2 are determined from the condition cii-Aj = lisbif A\ = x(7i/a)+y(7r/b) 
and A 2 = x(ir/a) — y{n/b) (fig. [7b). 

Since graphene is a covalent solid, a good approximation of its electronic band 
structure can be found using a tight-binding description |43] with two basis functions 
per unit cell (one p z orbital per carbon atom). This is justified in materials with sp 2 
bonding, in which the a bonds, localized in the graphene plane, are decoupled from the 
delocalized p z orbitals. Generalization from the unit cell to the entire solid is achieved 
by using Bloch function as ansatz. The band structure is described by: 

E = ±tyj 1 + 4 cos(k y b) cos(k x a) + 4 cos 2 (k y b) (8) 

The two bands (fig. [7b), a consequence of having two basis functions per unit cell, are 
symmetric about E = 0. In neutral graphene the Fermi level lies at E — 0, and all the 
states with E < are filled while those with E > are empty. In reciprocal space, 
the location of the charge-neutrality point (E = 0) is found by setting Eq. 0=0; with 
k x a = 0, the condition E = is satisfied at the six corners of the first Brillouin zone, 
shown in fig. UJc). The 6 corners fall in two groups of three equivalent points, differing 
only by a reciprocal lattice vector. These two groups represent then two non-equivalent 
points, named for crystallographic convention K and K'. 

Graphene is thus a zero band gap semi-metal, with the Fermi level located at the 
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intersections between the valence and conduction (ir and ir*) bands, which are located 
at the K and K' points in k-space. 

5.2. Low- energy spectrum 

The results from the tight-binding model are particularly interesting close to the K and 
K' points in reciprocal space, where the energy dispersion E oc k, in contrast to E oc k 2 
as in conventional electrical semiconductors. Around these K, K' points, the energy 
dispersion of carriers is similar to that of ultra-relativistic particles with zero rest mass 
mo, E(k) = \jm\d 1 + c 2 h 2 k 2 = chk; under these conditions the Schrodinger equation 
for Bloch electrons reduces to the 2D Dirac equation [HI H5]. The Dirac equation is 
thus used to describe the behaviour of carriers in graphene, which mimic massless Dirac 
Fermions around the K, K' (Dirac) points. 

By taking the K(K') point as reference and defining a vector momentum relative to 
this point as q = k — K (q' = k — K'), the eigenf unctions in momentum representation 
for the pseudo-relativistic carriers at the K and K' points are 

^,M = ^{e~ lSk/2 ±e^ 2 ) (9) 

^ Kl (q') = ^(e^' 2 ±e-^ 2 ) (10) 

with the +/— signs corresponding to energies of the ir* and ir bands. 

The wave functions at each Dirac point are thus described by a two component 
spinor, a linear combination of excitations arising from the A and B sublattice (ir and 
ir* bands). The pseudo-spin, quantified by the angle 9k, is in fact related to the vector 
momentum of the carriers and not to the real spin of electrons: 

9 k = arctan ( — J (11) 
\%) 

Since 9k is a function of momentum, it follows that the wavefunctions at K and K' are 
related by time- reversal symmetry (t — > —t). 

Another consequence of the pseudo-spin being associated to the vector momentum 
q, is that when carriers move in a full circle the wavefunction changes sign (or 
equivalently, it accumulates an effective phase e ±i7r ). Taking as an example one of the 
components of ^ +:K (k), a change of 360° = 2ir results in e -*( 6, ^+ 27r )/ 2 = e -w e -ie k /2 _ 
_ e -«ffe/2_ Thig j s a i so called geometric or Berry phase, and it has interesting implications 
whenever electrons in graphene move along closed trajectories, as for cyclotron motion at 
high magnetic fields or self-crossing paths in quantum diffusive regime at low magnetic 
fields. 

5.3. Half-Integer quantum Hall effect 

Novel and unique magnetotransport features are observed in graphene as a consequence 
of the facts that carriers: 1) can be modelled as having no mass, and 2) accumulate a 
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Berry phase of it when completing a full turn around the K, K' points. To start with, 
the expression for the cyclotron frequency for Dirac Fermions is modified from that in 
conventional 2D systems pQ: 



Additionally, instead of using the Schrodinger equation, the LL spectrum for chiral 
carriers in graphene is found by solving the Dirac equation in the presence of electric 
and magnetic field. The LL spectrum is given by [TJ H6j H7] : 

E LL „ Gr = ±%uj c Vn = v F V2keBN (13) 

with N an integer number including zero. The main differences with conventional 2D 
systems are: 

• The energy spacing of LL in graphene depends on the magnetic field as AE^l oc 
y/B, instead of AEn oc B as in conventional 2D systems. 

• Each LL in graphene can take twice as many electrons as LL do in conventional 2D 
systems. This four- fold degeneracy is due to spin- up/spin-down (as in conventional 
2D systems) , and valley degeneracy, K and K' (particular to monolayer graphene) . 

• The special LL in graphene at E = is shared equally by electrons and holes and. 
As a consequence it contains half as many states as the rest do. 

This anomalous Landau Level spectrum results in a peculiar sequence of plateaux 
in p xy in the quantum Hall regime, the so-called half-integer quantum Hall effect (fig. 

E}. 

P*y = ^( N + 1/2)" 1 = ^(4iV + 2)" 1 (14) 

with N an integer number including zero. This can be understood by recalling 

the case of conventional 2D systems, in which each (spin degenerate) Landau Level 

can allocate 2eB/h electrons (per unit area) (Eq. [5]). In addition to spin degeneracy, 

for graphene we need to take into account an additional double (valleys K and K') 

degeneracy. Thus, if we take into count only electrons (Fermi level Ep > 0), the electron 

density ng corresponding to N filled LL in graphene is: 

2eB, 2eB, 2eB , 4eB . 

n Gr = N—\ K + N—\ K . + -t-\e=o = -r-{N + 1/2) (15) 
n n n a 

By substituting the discretised density of states n<> into the B-dependence of 

the transversal resistivity p xy in conventional Hall effect, we arrive at the sequence 

of plateaux for the half-integer Quantum Hall effect: 

B B h 

Pxy en Gr e[4eB(N+ 1/2) /h] 4e 2 (N + l/2) 1 ' 

The half-integer quantum Hall effect is the fingerprint of monolayer graphene, and 

can be used experimentally to prove that electronic transport occurs through a single 

graphene layer. 
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Figure 8. (a) The sequence of Landau levels in graphene is unique; the energy spacing 
depends on the magnetic field as AEll oc v^B, instead of AEll oc B as in conventional 
2D systems and there exists a LL at E=0, shared equally by electrons and holes, (b) 
Comparison of quantum Hall effect (QHE, dashed lines) in conventional 2D systems 
and the half- integer quantum Hall effect (HiQHE, solid lines). For graphene, only 
v = 2 and v = 6 are shown for clarity. 



6. Quantum Hall effect in exfoliated graphene 

The accuracy of the resistance quantisation in graphene was first tested by measuring its 
value in terms of a room temperature 100 Q resistor with a CCC bridge described in the 
section [3TTT . The 100 Q resistor was measured prior and after the graphene measurements 
in terms of a standard well-characterised GaAs heterostructure device [12]. In such an 
indirect comparison the accuracy is ulimately limited by the short-term stability of the 
room-temperature resistor to a few parts in 10 9 . However, the first measurements by 
Giesbers et al. on a small exfoliated sample were at a much lower level of 15 part in 
10 6 . This was attributed in part to the small breakdown current the small device could 
sustain (2.5 pA for a width of 1 fim) which results in a small signal to noise ratio. One 
should note however that this small breakdown current still equates to a current density 
of 3.5 A/m which is equal or better than the best GaAs devices. A second limitation 
was the relatively high contact resistances (in the range of 1 kf2) which increases the 
noise in the measurements and leads to local heating [3]. 

The LNE group managed to improve on these initial results by measuring both 
exfoliated monolayer and bilayer graphene samples with significantly lower contact 
resistances (between 10-500 Q) [ISj. They made detailed measurements of both the Hall 
and longitudinal resistances and used the empirical relation, Ap xy = kp xx , to extrapolate 
t° Pxx — 0, the dissipationless limit for which perfect quantisation is expected. Both 
mono and bilayer graphene devices were found to be quantised to within 5 parts in 10 7 . 
This accuracy is again limited by the small breakdown current the exfoliated device 
could sustain (about 1 /j,A for a width of 2 fim). The LNE group determined that their 
devices suffer from a large degree of charged impurity scattering which could limit the 
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breakdown current and results in considerable charge density fluctuactions. 

Recently the PTB group managed to produce very large (30 /im by 150 /im) 
exfoliated graphene samples on top of a GaAs/AlAs heterostructure [19]. GaAs was 
chosen as a substrate because its surface is much smoother than the more commonly 
used SiC>2 surface and its higher dielectric constant should improve the electrical 
screening of substrate defects. It was also found that the GaAs surface results in 
much larger exfoliated flakes which, they speculate, might be as a result of the stronger 
hydrophilic character of GaAs. The heterostructure also contained a strongly doped 
GaAs layer which serves as a back gate to control the carrier density. The device had 
a breakdown current of ~ 15 fiA (i.e. a breakdown current density of 0.5 A/m) and 
all contact resistances were below 10 Q. The reported uncertainty for their precision 
measurements [50] was about 6 parts in 10 9 for zero back gate voltage which is on a par 
with the early precision measurements in epitaxial graphene of similar size [15]. The 
measurements were done at zero gate voltage because of a leakage between the graphene 
and back gate layer [19]. This leakage implies that a parallel conduction path exists to 
the graphene layer and this could lead to a systematic deviation from the true value of 
i?K which would need to be eliminated in future. 

It is very encouraging that accurate QHE results can be obtained with exfoliated 
graphene. However, the breakdown current and hence the achievable accuracy are at 
least an order of magnitude lower than those obtained in large area epitaxial graphene 
which is the topic of the next section. 

7. Graphene engineering for resistance metrology 

7.1. Growth and Characteristics of as- grown epitaxial graphene 

The highest quality graphene is produced by the well-known mechanical exfoliation 
from graphite (more commonly known as sticky-tape method) (4Tj. Although perfect 
for scientific research, their small size limits the applicability in high-precision metrology 
because of the need for high current densities and small contact resistances. Large-area 
graphene can be produced by two distinct methods, chemical vapour deposition (CVD) 
and high-temperature sublimation of silicon carbide. 

With the CVD technique very large areas of good quality graphene can be produced 
by decomposing and graphitising organic material on a hot metal surface. The metal 
acts as a catalyst for the reaction and typically copper is used as a substrate and methane 
as a carbon-containing gas. After growth, for most applications, the graphene needs to 
be transfered to a non-conducting substrate, a process which degrades the quality of 
the graphene by introducing doping and defects. Nevertheless, the Purdue group have 
demonstrated the QHE in a 7 mm by 7 mm graphene layer grown by CVD on copper, 
and then transferred to a Si02/Si substrate [5l] . Low fabrication costs are one of the 
main advantages of the CVD technique, while its challenges are the reliable control of the 
nucleation sites and the development of reliable, scalable and non-destructive graphene 
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transfer methods. Since the whole field is driven by the need for an indium-tin-oxide 
replacement by highly conductive graphene [52], the mentioned challenges must be met 
soon. However, in view of metrological applications, the low cost of graphene samples is 
of less importance, since it is negligible on comparison with the measurement equipment. 

Growth of graphene on SiC occurs when the crystal is annealed at high temperatures 
(< 1300 °C): the silicon sublimes and the carbon-rich surface re-crystalises to form 
graphene [53]. The quality and number of layers of graphene, critically depends on the 
precise growth conditions, crystal type, face and orientation. Compared to CVD, the 
electronic properties of the as-grown large area monolayer graphene produced by this 
method are often superior in terms of mobility and doping. The most commonly used 
polytypes of SiC are 4H-SiC and 6H-SiC, where H stands for hexagonal and the number 
refers to the number of basel planes (Si-C bilayers) in the unit cell. Both 4H-SiC and 
6H-SiC crystals have two different faces, either terminated by carbon or silicon atoms. 
The reaction kinetics of silicon desorption is different for the carbon and silicon faces, 
it is much faster on the C-face compared to the Si-face. 

Growth on the C-face leads to many graphene layers stacked on top of each other. 
These graphene layers are rotated with respect to each other and as result electrically 
decoupled to some extend and is therefore different from ordinary graphite [51]. Very 
high mobilities (of order 100000 cm 2 V _1 s~ 1 ) have been achieved in multilayer graphene 
but the QHE has never been observed in this system [55]. Recently the Georgia Tech 
group succeeded in growing a single monolayer graphene on the C-face of 4H hexagonal 
SiC [H] with mobilities in the range of 5000 — 20000 cm 2 V~ 1 s~ 1 and hole carrier density 
in the range of (1 — 2) x 10 12 cm -2 . A clear half integer QHE effect was observed in 
this system, although the metrologically relevant parameters were not explored in this 
work. 

Growth on the Si-face is much slower and allows better control of the thickness 
through the growth parameters; fabrication of monolayer graphene is possible on this 
face. A number of groups have succeeded almost simultaneously in observing the QHE 
by growing monolayer graphene on the Si-face in a high vacuum [13] and in an argon 
atmosphere [151 EH D2]- Typically mobilities in the range 500 — 5000 cm 2 V~ 1 s~ 1 are 
obtained and electron carrier densities in the range (0.5 — 10) x 10 12 cm -2 . Jobst et 
al. [IE] have made a careful study of the spread in device parameters (/i and n s ) by 
analysing data from over 50 devices. They found that for their growth temperature of 
around 1600°C and the 6H-polytype, the average mobility is 900 cm 2 V~ 1 s~ 1 and the 
carrier density is 1 x 10 13 cm -2 . This high carrier density makes it impossible to directly 
use these devices for quantum Hall effect resistance metrology (the v = 2 magnetic field 
would be larger than 100 Tesla). 

The SiC crystals are nominally cut at a 0° angle but in reality a slight miss-cut 
will lead to a number of atomically sharp step edges separated by flat terraces across 
the crystal surface. The Erlangen group investigated the effect of these step edges by 
placing small Hall bars either entirely on a terrace or deliberately across a number of 
step-edges. Surprisingly, there was no noticeable effect on the mobility [16]; theoretical 
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analysis of the role of steps on SiC on the charge carrier scattering in graphene, together 
with the review of experimental data is given in [56]. Growth parameters and degree of 
miss-cut determine the substrate morphology and the geometry of bilayer islands, which 
are nucleated along the substrate steps [57] and affect mobility 

Virojanadara et al. [58] found that the use of an inert gas during growth allows 
for an increase in the annealing temperature (~2000°C) by reducing the sublimation 
rate of silicon carbide. The advantage of the higher growth temperature is that it leads 
to atomically uniform graphene with larger graphene domains (~ 50 /mi). Tzalenchuk 
et al. [15] investigated the properties of large-area epitaxial graphene devices grown 
on the 4H-polytype using this technique. Twenty Hall bar devices of different sizes, 
from 160 /mi x 35 /mi down to 11.6 /im x 2 /mi were produced on each 0.5 cm 2 
wafer using standard electron beam lithography and oxygen plasma etching (Figured]). 
Atomic force microscopy (AFM) images revealed that the graphene layer covers the 
substrate steps like a carpet, preserving its structural integrity. Contacts to graphene 
were produced by straightforward deposition of 3 nm of Ti and 100 nm of Au through 
a lithographically defined mask followed by lift-off, with the typical area of graphene- 
metal interface of 10 4 /xm 2 for each contact. The manufactured material was n-doped, 
with the measured electron concentration in the range of (5 — 8) x 10 11 cm -2 , mobility 
about 2400 cm 2 V~ 1 s~ 1 at room temperature and between 4000 and 7500 cm 2 V~ 1 s~ 1 
at 4.2 K, almost independent of device dimensions and orientation with respect to the 
substrate terraces similar to results obtained in Ref. |16j . 

Figure |9] shows the longitudinal (dissipative) resistance, R xx and transverse (Hall) 
resistance, R xy of a 2 /mi wide Hall bar at 4.2 K and magnetic field up to 14 T. At high 
magnetic field we clearly identified two QHE plateaux, at R^J = R^/2 (N = 0) and 
R^j = Rk/6 {N = 1) corresponding to the filling factors v — 2 and v — 6, respectively. 
In graphene v = 2 corresponds to the fully occupied zero-energy Landau level (N = 0) 
characterised by the largest separation vyZheB (where v ~ 10 8 cms -1 is the Fermi 
velocity in graphene) from other Landau levels in the spectrum and hence the Hall 
resistance quantisation is particularly robust. This plateau appeared in the field range 
9-12 T, depending on the carrier concentration (which was beyond our control at this 
time) and was accompanied by a vanishing R xx . The N = 1 plateau, at v = 6, was not 
so flat, and R xx developed only a weak minimum. There was also a trace of a structure 
corresponding to v — 10. The observed sequence of Hall plateaux confirmed that the 
studied material was indeed monolayer graphene. At low magnetic fields Shubnikov-de 
Haas oscillations were observed as well as a weak localization peak characteristic of the 
phase coherence of electrons in a disordered system, analyzed comprehensively in [59]). 

The magneto-transport measurements on a much larger, 160 /im x 35 /im Hall bar 
device are also presented in fig. [9j A substantial positive magnetoresistance at low fields, 
which was absent in the smaller sample, indicated that the carrier concentration varied 
along the larger sample. Because of that, the v = 6 feature in R xx in the larger sample 
was less prominent. Nevertheless, despite the inhomogeneity of the carrier density, 
the Hall resistance plateau at R^j = -Rk/2 (n = 0) was accompanied by vanishing 
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Figure 9. Transverse (a) and longitudinal (b) resistance of a small and large device 
at T = 4.2 K with 1 current. (Adapted from Ref. [15] ) 

longitudinal resistance R xx . Importantly, the large-area device had a low contact 
resistance, R c ~ 1.5 Q, to the graphene layer and, as compared to smaller devices, 
could sustain a much higher current before QHE breaks down. Since larger breakdown 
current affords higher precision measurements in the QHE regime, the 160 fim x 35 fim 
Hall bar device has been chosen for the metrological measurements at the field of 14 T 
(limited by the size of our superconducting magnet, but still too far from the B = 17. 5T 
where the filling factor would be exactly v = 2 for this sample). 

The accuracy of the Hall resistance quantisation in graphene was established 
in measurements traceable to the GaAs quantum Hall resistance standard using a 
calibrated 100 Q resistor (see fig. [TO]). The optimal conditions at 300 mK were obtained 
for a source-drain current of 11.6 fiA, 15% below the breakdown current established in 
the measurements of R^. The quantisation accuracy +0.4 ±3 parts in 10 9 inferred from 
our measurements was a four orders of magnitude improvement on the previous best 
result in exfoliated graphene [12J. Graphene was still accurately quantised at 4.2 K, 
however, at this temperature the measurement current had to be reduced to 2.3 fiA, 
which increased the uncertainty of the data accumulated over a comparable time interval. 
The precision measurement presented in Ref. [15] readily puts epitaxial graphene 
quantum Hall devices in the same league as their semiconductor counterparts. Note 
that it was obtained on a sample, although large by graphene standards, substantially 
smaller than the semiconductor devices used for calibration and without optimisation. 
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Figure 10. Deviation between i?# and Rk in parts in 10 9 for a large graphcne sample. 
The accuracy is dominated by the data point at 11.6 fiA collected over 11 hours. The 
data points for higher currents are above the breakdown current and the device is no 
longer quantised. (Adapted from Ref. [15] ) 



Ideally, precision measurements are made at exactly integer filling factor because the 
breakdown current peaks at this value [3]. It is clear from the results presented in fig. [9] 
that this condition is not met for practical magnetic fields below 15 T. Further progress 
demands the ability to engineer the carrier density well below the 1 x 10 12 cm _2 -level 
and control the uniformity of charge distribution in epitaxial graphene (discussed in 
next section). 

Importantly, all twenty Hall bar devices fabricated on the chip shown in fig. Ufa) 
demonstrated quantum Hall effect, despite the fact that the steps on the substrate 
crisscrossed the devices. This tells not only that the graphene is continuous on 
these steps, but also that the bilayer patches, which are often nucleating along the 
substrate steps, did not cross even the narrowest 1 /im Hall bar. Bilayer patches 
crossing the sample cause scattering of the edge states and destroy the resistance 
quantization [60j EE]- Recent AFM studies show wide variations of the bilayer patch 
sizes and shapes, starting from submicrometer islands to a few micrometer long 
stripes [57J [62] • While the growth of single-layer graphene on the Si side of SiC is 
under control in many laboratories, the control of bilayer morphology is a more elusive 
task, and more work is required here. 

7.2. Carrier density control in epitaxial graphene 

There is a serious opportunity to use graphene in metrological applications. To reach 
the ideal v = nh/(eB) = 2 plateau, corresponding to the topologically protected 
N = Landau level at a feasible magnetic field, a low carrier density n in graphene is 
required. Therefore control of the carrier density in graphene is of utmost importance 
for metrology. 

Carrier density control in exfoliated graphene transferred onto a SiC^/Si substrate 
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Figure 11. Charge transfer between SiC and bare/gated graphene. (Adapted from 
Ref. [H]) 

is relatively easy. By applying a gate voltage to the doped silicon back gate, the carrier 
density can be changed by several orders of magnitude from electron to hole carriers 
right through the Dirac point |41j . A similar technique can be employed for graphene 
grown by CVD after transfer onto a an oxidised silicon wafer [51] . For epitaxial graphene 
grown on a highly insulating substrate this method is not possible and other top and 
bottom gating techniques have been developed to control the charge carrier density. 

As mentioned in the previous section, epitaxial graphene grown on Si-face of SiC 
is always strongly n-type doped with carrier densities typically in the range of 10 12 — 
10 13 cm -2 . The doping of the graphene is caused by the so-called "dead layer" or "buffer 
layer" of carbon atoms in between the SiC substrate and graphene (63J EH E51 EHJ EH EE] ■ 
This layer is non-conducting and characterised by a 6-\/3 x 6v3 supercell of the 
reconstructed surface of sublimated SiC. Missing or substituted carbon atoms in various 
positions of such a huge supercell in the dead layer create localised surface states with a 
broad distribution of energies within the bandgap of SiC. Kopylov et al. [19] developed 
a theoretical model of the charge transfer in epitaxial graphene to estimate the as-grown 
carrier density and degree of variation that can be induced by a top gate. 

The charge transfer results from the work function difference, A, between graphene 
and the combination of buffer layer and SiC substrate (fig. ITTT) . Assigning a density of 
states 7 to the donors in the buffer layer and pi to the donors in the bulk (p is volume 
density of donors and I the depletion length), the charge balance equation can be written 
as, 

7 [A — Aitde 2 n — e F ] + pi = n. 

(17) 

Here, ep = vhy/im, is the Fermi energy, and 47ie 2 dn results from the geometrical 
capacitance of the graphene layer. Typical values are A = 0.40 eV,d = 0.3 nm and 
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Figure 12. In direct chemical gating electrons in graphene can scatter on the ions 
in the adjacent active layer. When graphene is encapsulated in a polymer bilayer the 
active polymer (ZEP520A) is separated from graphene by a neutral spacer. (Adapted 
from Ref. HQ) 



7 = 5 x 10 12 eV _1 cm -2 . For an electrostatically gated device the charge is distributed 
between graphene and the gate, so we need to substitute n by (n + n g ) in Eq. [I7J 
n g = CV g /e is determined by the gate-to-graphene capacitance, C, and voltage, V g . 
Solving this equation [19] results in n ~ 10 13 cm -2 for zero gate voltage, which is close 
to the typical density of graphene doping observed in many recent studies of epitaxial 
graphene (57J Mi EH] . 

The above consideration suggest that the reduction of donors on or just under the 
SiC surface is crucial for the use of epitaxial graphene on SiC in gated devices. We note 
that for graphene samples grown at 2000°C [15] (a substantially higher temperature 
than in other published reports), reliably produced densities of n ~ 10 12 cm -2 much 
lower than 10 13 . Tzalenchuk et al. [10] speculate that high-temperature annealing of SiC 
compensates the surface donor states, possibly through surface segregation of impurities 
(e.g. B or N) abundant in bulk SiC crystals, or merely reduction of the number of surface 
defects. 

Top gate technology has been demonstrated for epitaxial graphene by a number of 
groups. Tanabe et al. [UJ produced a top gate stack of hydrogen silsequioxane (HSQ) 
and Si02 on top of a graphene Hall bar. The HSQ layer acted as an insulating layer 
and also protected the graphene from damage during processing. The metallic gate was 
formed by evaporating a layer of Cr/Au on top of the Si02- The mobility was unaffected 
by this processing and the carrier density could be changed from 10 12 cm -2 electrons 
through the Dirac point to 10 11 cm -2 holes. Shen et al. [13] used an 30 nm AI2O3 gate 
dielectric deposited via atomic layer deposition on a 1 nm oxidised Al seeding layer, to 
create a gate stack. In this device a modest 50% change in electron carrier density could 
be achieved. 

Back gate technology on Si-face epitaxial graphene devices was first demonstrated 
by the Erlangen group [70] , In this method a conducting layer is created some 700 nm 
below the surface of the SiC crystal through implantation of nitrogen ions prior to 
graphene growth. At low implantation dose and low temperatures this layer acts as back 
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Figure 13. Photochemical gating of graphcnc. (Adapted from Ref. |21j ) 

gate and the carrier density can be changed through the Dirac point if the graphene 
layer is decoupled from the buffer layer through intercalation of hydrogen. Without 
the intercallation step the donor states in the buffer layer will pin the Fermi level and 
effectively screen the graphene layer making the back gate ineffective. Recently, using 
the nitrogen implantation technique, Jounault et al. [7T| demonstrated the QHE for 
both electrons and holes in a bottom-gated epitaxial graphene sample on the C-face of 
SiC. Because the graphene layers in C-face SiC are almost decoupled there is no need 
for the hydrogenation step. 

A disadvantage of the gating technique is that an external voltage source has to be 
permanently connected to the device to maintain a constant carrier density. A solution 
would be static gate in analogy to semiconductor programmable nonvolatile memory 
devices. These devices are essentially transistors with one extra floating, isolated gate 
sandwiched between the control gate and the semiconductor channel. Charge can be 
transferred to the floating gate by an electric pulse on the control gate and stored 
there isolated almost indefinitely, until intentionally leaked through the dielectric, e.g. 
activated by UV light. In other implementations of the nonvolatile memory devices UV 
light is used for writing in which case thermal activation can be used for erasing. 

Lara-Avila et al. [21] demonstrated nonvolatile control of the carrier density 
by placing a polymer heterostructure of PMMA/MMA followed by another layer of 
ZEP520A, chosen for its ability to provide potent acceptors under deep UV light, on top 
of an epitaxial graphene device (both are commonly used clean-room photolithography 
resists). See the layout of the heterostructure and the chemical formulae of the polymers 
in fig. [121 the ZEP520A is the floating gate which is isolated from the graphene by the 
PMMA/MMA layer which acts as the spacer. In a sample with the initial carrier 
density n ~ 1.1 x 10 12 cm -2 , subsequent exposures to UV at 248 nm wavelength up 
to the dose 330 mJcm -2 decreased the low-temperature electron density 50 times down 
to 2 x 10 10 cm -2 (fig. [T3|) . This resulted in a fivefold increase in carrier mobility up to 
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Figure 14. Transverse (p X y) and longitudinal (yO xx ) resistivity measurement. The 
horizontal lines indicate the exact quantum Hall resistivity values for filling factors 
v = ±2 and ±6. (Adapted from Ref. [20]) 

16000 cm 2 V _1 s _1 at liquid helium temperatures and a tenfold increase in the resistivity 
of graphene. The irradiated devices remained latched in their high-resistivity state over 
many months. The "on/off" ratio of 10 for the resistivity of the photochemically-gated 
devices is similar to the best large-area single-layer graphene transistors demonstrated 
to date [72]. Very significantly, annealing the samples at 170°C, just above the glass 
transition temperature of the polymers, reversed the effects of light and returned the 
graphene charge carrier density to its value prior to UV exposure. 

Figure [H] shows the measurement of the longitudinal and transverse resistance as a 
function of magnetic field at low temperatures after encapsulation and UV illumination 
for the large Hall bar device (160 /xm x 35 /xm) and should be compared with fig. [9] It 
is clear that covering the device with the polymer sandwich has dramatically improved 
the homogeneity of the device besides reducing the carrier density. 

An alternative to photochemical gating is direct chemical gating as demonstrated 
by Jobst et al. [IB] . In this case the graphene is covered with tetrafluoro- 
tetracyanoquinodimethane (F4-TCNQ) molecules via thermal evaporation. Applying 
a monolayer of these acceptor-molecules reduced the carrier density down to 5.4 x 
10 11 cm -2 with an excellent mobility of 29000 cm 2 V s _1 . Although a good QHE was 
observed, the downside of this form of chemical gating is that it is not very stable in 
time, probably due to the absorption of water by the chemical layer. Recently, oxygen 
adsorption was used to reduce the carrier concentration in epitaxial graphene [73] and 
a quantum Hall effect has been observed at fields of 10 — 12 T, yet, the issues with the 
sample stability have still to be solved. 
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Figure 15. Schematic band-structure for graphene on SiC in zero field (a) and in 
quantising fields (b). (Adapted from Ref. [20] ) 

8. Magnetic field dependent charge transfer 

8.1. Filling factor pinning 

Figure [TTT b) shows magnetotransport data of a device in which the carrier density was 
reduced further down to 4.6 x 10 11 cm" 2 by photochemical gating [20]. An extremely 
wide v — 2 quantum Hall plateau as can be seen and the v = 6 plateau has all but 
disappeared in this case. A similar effect can also be seen in the low carrier density 
data of Ref. [16] . Janssen et al. [20] developed a model to explain this behaviour which 
is based on the pinning of the v = 2 quantum Hall state driven by a magnetic field 
dependent charge transfer mechanism. They argue that, specific to graphene on SiC, 
the pinning of the v = 4N + 2 filling factors is determined by the dominance of the 
quantum capacitance [73], over the classical capacitance, in the charge transfer between 
graphene and surface-donor states of SiC/G. 

The quantum capacitance of a two-dimensional electron system is the result 
of a low compressibility of the electron liquid determined by the peaks in j e . 
For electrons in high-mobility GaAs/AlGaAs heterostructures in magnetic field, the 
quantum capacitance manifests itself in weak magneto-oscillations of the electron 
density [751 EH] due to the suppressed density of states inside the inter-Landau level gaps. 
For epitaxial graphene on SiC, due to the short distance, d ~ 0.3 — 0.4 nm, between the 
buffer layer hosting the donors and graphene, the effect of quantum capacitance is much 
stronger, and the oscillations of electron density take the form of the robust pinning of 
the electron filling factor. Basically, the graphene layer and buffer layer exchange charge 
as the density of states oscillates with magnetic field. 

The charge transfer in SiC/G is illustrated in the sketches in fig. [151 for i? = 
(a) and quantising magnetic fields (b). In the latter case the Dirac cone transforms 
into discreet LL's which are depleted for successively larger magnetic fields. Again the 
transfer can be described using the charge balance Eq. [TTJ discussed earlier [T9| and the 
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Figure 16. Graphical solutions for the carrier density as a function of magnetic flux 
density, n s (B), of the charge-transfer model given by Eq. [17] (black line) together 
with lines of constant filling factor (red lines) and n s (B,N) (green lines) for n g = 
5.4 x 10 11 cm" 2 (a) and n g = 8.1 x 10 11 cm" 2 (b). (Adapted from Rcf. [20]) 



magnetic field dependent solution is shown in fig. [16] Interestingly, in regime II where 
the chemical potential in the system lies inside the gap between N = and N = 1 
LL, the carrier density increases linearly with the magnetic field, n s = 2eB/h, due to 
the charge transfer from SiC surface. The final carrier concentration of graphene is 
determined mostly by the density of donor states 7 in the buffer layer and can be nearly 
30% higher than the zero- field density |20j . 

Accurate quantum Hall resistance measurements require that the longitudinal 
voltage remains zero to ensure the device is in the non-dissipative state, which can 
be violated by the breakdown of the QHE at high source drain current levels. Fig. UTT a) 
shows the determination of the breakdown current I c for different values of the magnetic 
field along the v = 2 plateau. Here I c is defined as the source-drain current, I s d, at which 
V xx >10 nV [fig. Mb)]- 

From fig. [TTT b) we see that the breakdown current continues to increase in the v — 2 
quantum Hall state reaching m 500 /xA at the maximum field of 14 T. This behaviour 
is very different from that observed in conventional semiconductor systems where the 
breakdown current peaks at the exact integer filling factor [3] indicated by the dashed 
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Figure 17. (a) Measurement of V xx as a function of source-drain current at different 
values of magnetic field ranging from 7 T to (left hand curve) to 14 T (right curve) in 
steps of 1 T. (b) Transverse (p xy ) and longitudinal (p xx ) resistivity measurement at 
the reduced carrier density of4.6xlO n cm~ 2 measured at I s d = 1 /iA together with 
the measured break-down current, I c . Dashed blue line indicates position of the exact 
v = 2 filling factor for the low field carrier density. (Adapted from Ref. [32] ) 



blue line in fig. [TTT b). In graphene the filling factor is effectively pinned at v — 2 over 
a broad range of magnetic field resulting in a novel quantum Hall state which is ideally 
suited for high precision resistance metrology. The anomalous pinning is responsible for 
the extremely high breakdown current and wide operational parameter space of epitaxial 
graphene. 

When the Fermi level lies in between the N = and 1 Landau levels the activation 
energy for the dissipative transport hyjl/2v /\g is very large ~ 1000 K [77]. For such 
a high activation energy, the low-temperature dissipative transport is most likely to 
proceed through the variable range hopping between surface donors in SiC involving 
virtual occupancy of the Landau level states in graphene to which they are weakly 
coupled. Indeed, as shown in fig. ITSTa). the temperature dependence of the conductivity 
a xx measured obeys an exp(—^T*/T) dependence typical of the hopping mechanism. 
The T* values determined from the measurements at different magnetic fields are 
plotted in fig. fTBT b). The breakdown current rising with field to very large values 
[fig- flTT a)] corresponds to T* reaching extremely large values in excess of 10 4 K - at 
least an order of magnitude larger than that observed in GaAs [78] and more recently in 
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Figure 18. (a) Variable-range hopping analysis of the data in fig. fTTl b) for the 
magnetic field range from 5 to 14 T and temperature range from 30 to 100 K. The 
axis are rescaled to give a straight line fit for a xx (T) cx (1/T) exp (^— ^/T*/t) . (b) 
Hopping temperature T* as function of magnetic field. 



exfoliated graphene [791 ED] • Such an extremely large hopping temperature is difficult 
to explain within the existing theoretical models and will need further investigation to 
fully understand. 

8.2. Robustness of the quantum Hall state 

Robustness of the quantum Hall state is important for resistance metrology and implies 
that the quantisation is insensitive to variations in experimental parameters making 
it easy to realise the correct value of h/Ne 2 . Figure [19] demonstrates the robustness 
as a function of magnetic flux density for the v — 2 quantum Hall state [32]. Here 
Ap xy = 2[pg aphcnc ( J B) - pGaAs/AiGaAs^ = m5 T))/R K measures the difference between 
graphene and a fixed reference GaAs/AlGaAs device. The magnetic field for the 
GaAs/GaAlAs device is held constant at exactly v — 2 and the magnetic field for 
the graphene sample is varied. In a separate measurement the longitudinal resistivity 
was measured for the graphene device and plotted in the same graph. The graph shows 
that the quantum Hall effect is accurate over a range of at least 4 T and only limited 
by the maximum magnetic field. The inset is a plot of Ap xy as a function of p xx and 
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Figure 19. Measurement of Ap xy = 2[p^ a P hcno (B) - p° aAs/A1GaAs (B= 10.5 T)]/R K 
and p xx as a function of magnetic flux density on the graphene device. T = 0.3 K 
for the graphene device and 1.5 K for the GaAs/AlGaAs device. The measurement 
current was 60 /iA. ler-error bars for Ap xy are of the order of 5 ppb and not visible 
on the scale of this graph. The inset shows a plot of the same data but here —Ap xy is 
plotted as function of p xx . Blue line linear fit. (Adapted from Ref. [321) 



confirms the empirical relationship Ap xy = kp xx with k = 0.39 which has been observed 
for traditional semiconductor systems [3] (The numerical value of k is not universal and 
different for different contact pairs and different cool-downs of the device [32]). 

In figure [20] the same experiment is repeated but this time with temperature as 
the parameter at a magnetic field of 14 T. The graphene device shows good (ppb-level) 
quantisation up to at least 15 K, however it is worth noting that at the our maximum 
magnetic field we are probably not yet at the plateau center and so the maximum 
operational temperature could be higher. The inset shows the plots of Ap xy versus p xx 
for opposite magnetic field directions from which we determine that k = 0.16. The sign 
of k changes with the magnetic field direction is again in accordance with that observed 
for semiconductor systems. A k value of 0.16 implies that a relative error in R^=2 hcn ° < 
parts in 10 10 is achievable when p xx < 10 pfl. 

9. Energy loss and breakdown 

The large breakdown currents observed in graphene mean there is significant heating of 
the electron gas in the quantum Hall state (for 100 pA this is more than 100 pW in a 
35 pm wide channel). This in turn implies that the charge carriers in graphene must be 
very efficient in losing their energy to the lattice, a property which makes this material 
so attractive for modern electronics where the power densities are limiting performance. 
This property is somewhat counter intuitive given that the available phonon spectrum 
for scattering is very limited (up to room temperature polar-optical phonon scattering 
and piezoelectric acoustic phonon scattering are negligible). 

Baker et al. [EH [82] studied the energy loss rates of hot charge carriers in graphene 
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Figure 20. Measurement of Ap xy = 2[p^ a P hene (T) - p^ As/MGaAs (T= 1.5 K)]/R K 
and p xx as a function of temperature for the graphene device. B = 14 T for the 
graphene device and 10.5 T for the GaAs/AlGaAs device. The measurement current 
was 60 /iA. lcr-error bars for Ap xy are of the order of 5 ppb and not visible on the 
scale of this graph. The inset shows a plot of the same data but here Ap xy is plotted 
as function of p xx for opposite magnetic held directions. Red and green lines linear 
fits. (Adapted from Ref. (32]) 



produced by three different methods, exfoliated, epitaxial and CVD-grown samples. 
By measuring the amplitude of the Shubnikov-de Haas oscillations as a function of 
temperature and current, the carrier temperature can be determined as a function of 
current. For epitaxial graphene in which not enough SdH oscillations are present, the 
amplitude of the zero field weak localisation peak was used to determine the carrier 
temperature. From these results the energy loss rates can be obtained as a function of 
carrier temperature which are shown in fig. |21~1 

The data in fig. [21] show a behaviour which is consistent with a ~ T 4 power law 
dependence. The solid line is based on a theoretical model by Kubakkadi for acoustic 
phonons, extrapolated from low temperatures [83]. Using this theoretical model, the 
electron energy relaxation time, r e , can be obtained from the energy loss rate. Baker et 
al. typically find that r e in graphene is at least an order of magnitude smaller than that 
observed in GaAs over a wide temperature range. This is despite the fact that optic 
phonon emission is making a strong contribution to the GaAs energy loss rate above 
~ 50 K. 

This high loss rate together with the large cyclotron energy explains the high 
breakdown current in graphene. This can be seen from the most commonly used model 
to predict the breakdown of QHE, which is the bootstrap-type electron heating model 
of Komiyama and Kawaguchi [M] , which is based on the runaway heating which occurs 
when the quantum Hall effect begins to break down. Within this model, for graphene, 
the breakdown field E y is predicted to be 



E y 



1 Bhcu c 
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Figure 21. Carrier energy loss rate as a function of electron temperature for a 
representative variety of samples and from the two techniques used. The data follow 
a similar trend for all samples, with the energy loss rate scaling by approximately T4. 
An extrapolated T 4 dependence from low temperature calculations |83j is shown for 
a carrier density of 4.13 x 10 11 cm" 2 and a T 3 dependence is also shown. (Adapted 
from Rcf. [82]) 

We see that with typical parameters for graphene (r e = 4 ps and hu c = 126 meV) 
we can expect a breakdown current 10 times higher than for GaAs (r e = 100 ps and 
hu c = 17 meV) [81]. These results suggest that the intrinsic properties of graphene 
make it a much better choice for the realisation of a quantum resistance standard than 
traditional GaAs. 

10. Universality of the QHE in epitaxial graphene 

As discussed in section [3] in order to test the universality of the quantum Hall effect, 
two different QHE devices need to be set up at the same quantum Hall plateau so that 
a one-to-one comparison of resistance can be made (in principle one could also compare 
different index plateaux if different winding ratios are available on the CCC). In 2010 
NPL and BIPM collaborated in a joint experiment to compare a graphene device and a 
GaAs/AlGaAs device directly against each other by using the BIPM traveling quantum 
Hall system |22j . 

The graphene sample was mounted in a 14 T/300 mK cryostat and connected to 
the Slave side of the CCC bridge. Two GaAs/AlGaAs samples were mounted in the 
transportable 11.7 T/1.2 K BIPM cryostat and connected to the Master side. The 
two GaAs samples used were traditional GaAs/AlGaAs heterostructures obtained from 
the PTB (device 1) and LEP (device 2). The red triangles in fig. [22] are the results 
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Figure 22. Measurement of AQ a A s /AiGaAs- Graphene expressed in parts per billion 
for the direct comparison of /j^aAs/AiGaAs anc [ ^Graphene &g a f unc ti on f j gd "pjjg 
uncertainty bars represent the ±ler random error of the mean. Different symbols are 
explained in the text. The red line is the weighted mean of all the data points and the 
green lines signify ±lcr. (Adapted from Ref. \'22\ ) 



for graphene against device 1 for 4 different source-drain currents in the devices. The 
pink diamond is a measurement for graphene against device 2 at at I S( i = 50 fiA. Here 
A G aAs/AiGaAs-Graphenc = (^R GaAs / MGaAs _ #| raphene )/(# K /2) and each data point consists 
of an average of between 3 and 10 hours worth of data. The uncertainty increases for 
lower I sd because the signal-to-noise reduces for lower I sd . 

To check for errors due to non-zero p xx the measurements were repeated for non- 
opposite contacts (green dot and blue square). (Note that it is very difficult to measure 
p xx directly to the required level of precision.) Another test to check for small errors 
is to reverse the direction of magnetic field on the graphene sample. The results of 
this measurement is represented by the light blue hexagon. Finally, the devices were 
exchanged between the NPL and BIPM cryostats in order to eliminate small parasitic 
leakages (black square in fig. 122]) . 

The data in fig. [22] show no sign of any systematic errors in the measurement 
campaign and so all results can be combined to give a weighted mean of 
AcaAs/AiGaAs-Graphene = (-4.7 ± 8.6) x 10 -11 . The random noise of 8.6 parts in 
10 11 dominates over the other components, estimated to have a combined standard 
uncertainty of 1.6 parts in 10 11 |32j . 

Previously our knowledge of the universality of the QHE has been limited to the 
level of 2 or 3 x 10~ 10 for comparisons between GaAs and Si or between identical GaAs 
devices [SI [39] [33] [34] . However both GaAs and Si are traditional semiconductors with 
a parabolic bandstructure and governed by the same physics. Graphene is a semi metal 
with a linear bandstructure and is described by Dirac-type massless charge carriers 
and so universality in terms of material independence goes well beyond the comparison 
between two semiconductors. The result does directly support the Thouless-Laughlin 
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argument [6] that the Hall conductivity is a topological invariant and is a fundamental 
test of condensed matter theory. 

The value of A GaA s/AiGaA s -Gra P h e n e = (-4.7 ± 8.6) x KT 11 on the material 
independence is the strongest evidence yet of the hypothesis that the resistance is 
quantised in units of h/e 2 is correct. This result underpins the assumptions made 
in the Watt balance experiment and thereby supports the pending re-definition of the 
Si-units for kilogram and ampere in terms of h and e. 

11. The era of graphene-based metrology is only just beginning 

11.1. Graphene single- electron pumps for better quantum current standards 

Single electron pumps are set to revolutionize electrical metrology by enabling the 
ampere to be re-defined in terms of the elementary charge of an electron [85]. Pumps 
based on lithographically-fixed tunnel barriers in mesoscopic metallic systems [861 EZ] 
and normal/superconducting hybrid turnstiles [EH [89] can reach very small error 
rates, but only at MHz pumping speeds corresponding to small currents of the order 
1 pA. Tunable barrier pumps in semiconductor structures have been operated at GHz 
frequencies [90j[9T], but the theoretical treatment of the error rate is more complex and 
only approximate predictions are available [92]. 

Connolly et al. [93] have recently demonstrated a monolithic, fixed-barrier pump 
made entirely in graphene. Their adiabatic double quantum dot pump operates at 
frequencies up to 1.4 GHz with a predicted error rate as low as 0.01 parts per million 
at 90 MHz. Ten pumps operated in parallel would deliver 100 pA with metrological 
accuracy. The authors ascribed the high performance of the pump to a number of factors 
deriving from graphene's unique two-dimensional physical and electronic structure: 

Firstly, the presence of strong edge and potential disorder in lithographically 
defined graphene nanostructures leads to the formation of multiple quantum dots in the 
constrictions acting as tunnel barriers between the dots. Rather than this impeding 
pumping, the resulting capacitance and resistance of the random tunnel junctions 
between the quantum dots promotes a high intrinsic tunneling time while simultaneously 
suppressing co-tunneling events due to the large overall dissipation. 

Secondly, the large interdot capacitance and linear electronic dispersion in graphene 
leads to a large and occupancy-dependent single-particle energy spacing near the Dirac 
point, A(iV) = hvp/(d\/N) (N ^> 1 - the number of electrons on the dot at a 
given Fermi energy, d - the dot diameter), which suppresses photon assisted interdot 
transitions and protects adiabaticity of the pump when operated at high frequency. This 
permits the interdot capacitance to be optimized through increased dot size, thereby 
improving the tolerances for parallelization, without compromising on accuracy. 

Combined with the record-high accuracy of the quantum Hall effect and proximity 
induced Josephson junctions [M], accurate quantized current generation brings an all- 
graphene closure of the quantum metrological triangle within reach [95J . 
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11.2. Quantum Hall arrays for resistance metrology 

The value of the quantum Hall resistance at v = 2 of « 12.906 kSl is not very practical for 
everyday resistance metrology which is concentrated around decade values of resistance. 
A large ratio, typically more than 100 to 1, has to be bridged in order to calibrate a 
100 ft resistance which is a demanding measurement. A solution to this problem is the 
use of quantum Hall arrays with which a wide scale of resistance can be covered [96], [97]. 
The use of QHARS has so far been limited. This is primarily due to low tolerance 
of Hall quantisation in the semiconductor arrays to any spread of characteristics in 
individual devises leading to low overall breakdown currents and hence low signal-to- 
noise ratio. In addition, connecting several hundred Hall bars in a series or parallel 
network, which requires many bridging connections and all contact resistances to be low 
Ohmic, entails complex multilayer fabrication. In contrast, large area graphene offers a 
wide operational parameter space for the quantum Hall effect, high breakdown current 
and the ease with which low Ohmic contacts can be achieved. Quite remarkably, a series 
array can be produced based on the ambipolar nature of graphene without the need for 
complex multi-layer interconnects [98J. 

11.3. Constant light absorption in a broad spectral range for radiometry 

Graphene is an attractive material for use in optical detectors because it absorbs light 
from mid-infrared to ultraviolet wavelengths with nearly equal strength, « not per layer 
[99] . Because of its small heat capacity, very weak coupling to thermal baths and rapid 
thermalization |100[ [82] it is particularly well suited for use in cooled bolometers - 
devices that measure the power of incident electromagnetic radiation via the heating 
of material with a temperature-dependent electrical resistance. Jun Yan et al. [101] 
have demonstrated a graphene bolometer, which exhibits a noise-equivalent power 
33 fW Hz -1 / 2 at 5 K that is several times lower, and intrinsic speed > 1 GHz at 10 K 
that is three to five orders of magnitude higher than commercial silicon bolometers 
and superconducting transition-edge sensors at similar temperatures. Making use of 
wide-bandwidth noise thermometry, Fong et al. |102] were able to measure temperature 
oscillations with a period of 430 ps. Modern absolute radiometers |103j . devices that 
measure the radiant power of electromagnetic radiation, are used to realise the candela 
- the SI unit of luminous intensity. They are based on the substitution principle, 
whereby the power of the radiant heating of an absorbing material is compared to the 
electrical heating, and operate at cryogenic temperatures. Although the absorptivity can 
reach 99.998% it is never exactly known leading to uncertainties in the measurements. 
Currently primary standard cryogenic radiometers can measure the power of an intensity 
stabilised laser beam to about 0.01%. Combination of the accurately known (though 
relatively small) wide-band optical absoption in graphene and high sensitivity and speed 
of graphene bolometers, can lead to novel quantum-based, fast, fairly compact and 
simple radiometers, linking the optical scale to the fine structure constant. 
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11. 4. Exceptional mechanical properties for displacement, force and mass sensors on 
the nanoscale 

Although the most precise mass sensing experiments with a resolution of 1.7 yoctogram 
(1 yg = 10~ 24 g), which corresponds to the mass of one proton, have been performed 
using carbon nanotube NEMS |104j . replacing them with graphene |105j brings the 
advantage of a larger area for force application and better coupling to the displacement 
sensor, although compromising the resonant frequency and hence the quality factor. 
The fundamental problem of the resonance-based mass (or force) measurements is that 
one cannot tell the mass (force) distribution on the resonator from the frequency shifts 
(one cannot hear the shape of the drum [106] ) even when multi-mode measurements 
are carried out. In the case of mass sensing, this problem can potentially be solved 
by functionalising graphene or nanotube, so that additional mass (e.g., a gas molecule) 
is absorbed at an atomically defined position. For force sensing graphene NEMS can 
be driven by the local radiation pressure and read-out also optically. Having in mind 
the quantised optical absorption of graphene this may give a unique way of realising 
sensitive force measurements linked to fundamental constants. In addition, graphene 
strain gauges with either electrical or optical readouts look very promising. Graphene 
is the only crystal which can be stretched by 20% |107j . thus enhancing the working 
range of such sensors significantly. 

11.5. Extreme sensitivity of electron transport for bio-molecular and gas metrology 

Detection of individual events when a gas molecule attaches to or detaches from 
graphene's surface have been demonstrated |108j . The adsorbed molecules change 
the local carrier concentration in graphene one by one electron, which leads to step- 
like changes in resistance. The achieved sensitivity is due to the fact that graphene 
is an exceptionally low-noise material electronically. This opens up unprecedented 
applications requiring counting of species, such as the measurements of the amount of 
substance on the nanoscale. For analytical applications however these measurements 
should be made discriminatively between different atomic/molecular entities, which 
requires chemical functionalisation of the graphene sensor. 

We can confidently predict that other, even more exciting applications of graphene 
in metrology will be investigated. It can be argued |109j that the major advantage 
of graphene sensors is their multi-functionality. A single device can be used in 
multidimensional measurements (for example, strain, gas environment, pressure and 
magnetic field). Whether or not graphene will fulfil the expectations is yet to be seen. 
But if we can practically realise the potential of this unique two-dimensional material, 
than the impact on technology and metrology will no doubt be revolutionary. 
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12. Conclusions 

The journey from original discovery of the QHE in graphene to a quantum resistance 
standard which outperforms the established technology in many aspects has been 
remarkably short. For epitaxial graphene the robustness of the quantisation in terms 
of temperature, magnetic field and source-drain current is exceptional. The material 
is cheap and relatively easy to fabricate and process. It allows for the realisation of a 
quantum resistance standard with modest means, e.g. a small superconducting magnet 
and cryocooler. As such it will improve the proliferation of quantum standards and 
allow many smaller laboratories to realise their own resistance scale. One even could 
envisage university students being able to perform QHE experiments, much in the same 
way as the discovery of High-T c superconductors enabled table-top experiments with 
Josephson junctions in many science classes. 
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